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Abstract. We study the two-dimensional ±J Ising model, three-state Potts model 
and four-state Potts model, by the numerical transfer matrix method to investigate the 
behaviour of the sample-to-sample fluctuations of the internal energy on the Nishimori 
line. The result shows a maximum at the multicritical point in all the models we 
investigated. The large sample-to-sample fluctuations of the internal energy as well as 
the existence of a singularity in these fluctuations imply that the bond configuration 
(or, equivalently, the distribution of frustrated plaquettes) may be experiencing a non- 
trivial change of its behaviour at the multicritical point. This observation is consistent 
with the picture that the phase transition at the multicritical point is of geometry- 
induced nature. 



PACS numbers: 05.50.+q, 75.50.Lk 
1. Introduction 

The two-dimensional spins glasses, especially the ±J Ising model and the random 
g-state Potts model, have been studied both by theoretical and numerical methods. 
Some exact results have been derived Q, the phase diagram has been drawn [|], § and 
recently the exact location of the multicritical point has been predicted [[|. However 
the characteristics of the multicritical point are not well understood except that the 
critical exponents have different values than those along the ferro-para critical curve 
|| IJ. The properties of the multicritical point are of great interest because this point 
represents one of the simplest non-trivial renormalization group fixed points in random 
spin systems [[J and also because we may be able to understand the significance of the 
Nishimori line from the study of the multicritical point. 

In the present paper we investigate the sample-to-sample fluctuations of the internal 
energy near the multicritical point. The thermal average of the Hamiltonian, namely 
the internal energy, has different values from one bond configuration to another, given 
a fixed probability parameter (p for the ferromagnetic bonds in the ±J model, for 
example). Although the internal energy is self-averaging, the variance of its distribution 
is of order N, the system size. If the variance is anomalously large and singular at 
some point, it would mean that a non-trivial (and possibly drastic) change happens in 
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the bond configuration at that point because the value of the internal energy fluctuates 
strongly from sample to sample. This is exactly what we shall show in the present paper 
to exist at the multicritical point by observing the system along the Nishimori line. 

The Nishimori line plays an important role in the determination of the phase 
diagram. In particular it is believed that the multicritical point lies on the line We 
have another advantage to study the sample-to-sample fluctuations along the Nishimori 
line: These fluctuations are directly related to the average specific heat on the line. This 
fact reduces numerical efforts considerably because it is not necessary to generate very 
many samples to calculate the variance. Only the average of the specific heat needs to 
be computed. We shall make full use of this relation later. 

In the next section, we show that the sample-to-sample fluctuations of the internal 
energy are very simply related to the specific heat on the Nishimori line. This facilitates 
numerical calculations enormously as mentioned above. In section 3 we investigate 
numerically the behaviour of these fluctuations near the multicritical point on the 
Nishimori line for three models: the ±J Ising model, three-state Potts model and 
four-state Potts model. We will see that these fluctuations show a maximum at the 
multicritical point. The physical significance of the results is discussed in the final 
section. 



2. Average fluctuations of the internal energy 

Let us consider the random Z q model with gauge symmetry which includes the ±J 
Ising model and the random g-state Potts model. The general Hamiltonian is written 
as follows: 

H = - J £V(S i -S j + J ij ), (1) 

(ij) 

where Si = 0, 1, q — 1 is the g-state spin variable, Jy = 0, 1, q — 1 is the quenched 
random bond variable and V is an interaction function of period q. The summation is 
taken over the first-neighbouring sites of the square lattice. The variable follows the 
probability distribution 

P(J^P P ) = - ] • (2) 



9-1 

E 



c 

77=0 

The Hamiltonian (1) is invariant under the gauge transformation Si —>■ Si — o~i, Jij — > 
Jij + o~i — aj, where cr, = 0, 1, q — 1, but P(Jij] (3 P ) is not. Using the gauge theory [[|] 
some exact results can be derived such as the average internal energy on the so-called 
Nishimori line defined by /3 = 1/T = (3 pi 

[E] = [(H)] = -N B \ ln j v{v) . (3) 
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Here the square brackets denote the configurational average, the angular brackets are 
for the thermal average and N B is the number of bonds. If we denote 

77=0 

equation (3) is expressed as 

[E] = -N B {V}. (4) 

In the case of the ±J Ising model, q = 2, V(0) = +J and V(l) = —J, the average 
internal energy (3) reduces to 

[E] = -iV B Jtanh(/3J). 

In the case of the g-state Potts model with V(0) = +J and V(l) = V(2) = ... = 
V(q-1) = 0, 

\E\ = -N B J — — -rj. 

[ J l + (q- l)e"^ J 

We also find from the gauge theory that 

[(H 2 )] = N B {V 2 } + N B (N B - 1){V}\ (5) 
when p = P p , thus, with equation (4), 

Nb({V 2 } - {V} 2 ) = \(H 2 )} - [(H)] 2 . (6) 
Since the specific heat is defined by 

T 2 [C] = [(H 2 ) - (H)\ (7) 
we obtain, combining equations (6) and (7), 

Nb({V 2 } - {V} 2 ) - T 2 [C] = [E 2 ] - [E] 2 (> 0). (8) 

The right-hand side of equation (8) is the sample-to-sample fluctuations of the internal 
energy, [(AE) 2 ]. Note that equation (8) holds only on the Nishimori line f3 = f3 p . In the 
case of the ±J Ising model equation (8) reduces to the well-established inequality 

T 2 [C] < N B J 2 sech 2 {(3J). 

Equation (8) indicates that the sample-to-sample fluctuations of the internal energy 
can be evaluated from the average of the specific heat because the quantity {V 2 } — {V} 2 
is trivially calculated. This fact is a great advantage since it eliminates the necessity 
to evaluate the variance directly by using a large number of samples; the specific heat 
is a self-averaging quantity and therefore can be calculated with high precision from a 
small number of samples with sufficiently large system size. Away from the Nishimori 
line, we do not have such a simple relation because equations (4) and (5) do not hold. 

Since the actual distribution of bond values (and therefore of frustration) changes 
from sample to sample in finite-size systems with given /3 P , the internal energy takes 
various values from one set of bonds to another around the average value. The relation 
(8) helps us to estimate the variance of the internal energy [(AS) 2 ], and directly depends 
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Figure 1. Phase diagram of the ±J Ising model in two dimensions. There exist 
ferromagnetic (F) and paramagnetic (P) phases, and no spin glass phase exists. The 
Nishimori line is shown dashed. 

on the distribution of frustration, which has a geometrical nature: The sample-to- 
sample fluctuations of the internal energy should reflect the geometrical structure of 
the distribution of frustration, which is a quantity independent of thermal variables 
but depends only on the bond configuration. This argument suggests that we may 
be able to investigate the properties of the distribution of frustration by looking at 
the sample-to-sample fluctuations of the internal energy, or equivalently, the average 
specific heat on the Nishimori line. It is believed that the Nishimori line intersects 
the phase boundary between paramagnetic and ferromagnetic phases [Q (see figure []]). 
If the distribution of frustration goes under a sharp change of its behaviour, it would 
be reflected in anomalously large sample-to-sample fluctuations of physical quantities 
including the internal energy. Such an anomaly, if it exists on the right-hand side of 
equation (8), is related to the singularity of the specific heat at the multicritical point 
on the Nishimori line due to equation (8). Since the specific heat is singular at the 
multicritical point, we may conclude that the sample-to-sample fluctuations are also 
singular there. Hence, by calculating the average specific heat near the multicritical 
point, we are able to investigate the structure of the distribution of frustration, which is 
of purely geometrical nature. A similar argument has been used to relate the free energy 
and the entropy of frustration |]]. The next section gives us some numerical supports 
to confirm this idea. 

3. Numerical calculations 

We have tested the ideas of the previous section by numerically calculating the average 
specific heat on long strips by transfer matrix method. We have studied three models: 
the ±J Ising model, three-state Potts model and four-state Potts model. We have 
used long strips of length 5000 and of various widths with free boundary conditions 
and without external field. By the transfer matrix method we calculated the free 
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energy along the Nishimori line. The specific heat was obtained from numerical second 
derivative of the free energy by the inverse of the temperature. We averaged the specific 
heat over ten samples. 

For the ± J Ising model the widths of the lattice are 3 < L < 14. The lattice size (up 
to 5000 x 14) and number of samples (ten) are much smaller than the existing numerical 
transfer matrix studies [|], HJ. However, since our purpose is much more modest than 
these investigations (not trying to evaluate the critical exponents, for instance), our 
system size turns out to be sufficient as we shall show. Another reason for the small 
scale of our computation is that the use of the formula (8) greatly reduces the necessary 
number of samples as discussed in the previous section. 

We define the parameter p for the ± J Ising model by 

e 2/3 p J = P 

1 — p 

Thus p is the probability that Jy = i.e. V(Jij) = +J |j. Each bond follows the 
probability distribution 

P(J ij )=p8(J ij ) + {l-p)S(J ij -l), 

and we restrict our studies to 0.5 < p < 1. For the three-state Potts model, 3 < L < 9 
and for the four-state Potts model, 3 < L < 7. For the g-state Potts model, we define 

e^ J = --(g-l), 
p 

and the probability distribution (2) is then 

p(j {j ) = (i - (q - i)p)5{j i3 ) + - I), 

1=1 

with < p < l/(q - 1). 

For each point on the Nishimori line f3 = (3 P we extrapolate the specific heat to 
the limit L — > oo. We have found that the extrapolation by the L _1 -law works quite 
well in almost all cases (see figures |^, ^] and This is natural because we are mostly 
treating values of p away from the multicritical point. Exactly at the multicritical 
point, the extrapolation needs some care, in particular if one wishes to estimate the 
critical exponents Q. However, for our purpose to confirm the existence of a peak 
in the sample-to-sample fluctuations of the internal energy at the multicritical point, 
the method of extrapolation is rather irrelevant, which we confirmed by using other 
powers such as L _L5 in some cases. The error bars in the extrapolated values have been 
estimated by changing the range of extrapolation fitting, for example, from comparison 
of extrapolation using 3 < L < 14 with the one using 6 < L < 14. 

The resulting values of the energy fluctuations are plotted in figures |, | and |7| for 
the ±J Ising model, three-state Potts model and four-state Potts model, respectively. 
The position of the multicritical points are predicted in reference 0: for the ±J Ising 

| Notice that we use the Hamiltonian (1) where Jy takes the values and 1. This notation is completely 
equivalent to the usual Hamiltonian H = — JijSiSj where Si,Sj and Jij / J take the values ±1. 
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Figure 2. Extrapolation of T 2 C near the multicritical point for the ±J Ising model, 
p = 0.8900. 
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Figure 3. Extrapolation of T 2 C near the multicritical point for the three-state Potts 
model, p = 0.0795. 
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Figure 4. Extrapolation of T 2 C near the multicritical point for the four-state Potts 
model, p = 0.0615. 
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Figure 5. Sample-to-sample fluctuations of the internal energy for the ± J Ising model. 



model p c = 0.889972, for the three-state Potts model p c = 0.0797308 and for the four- 
state Potts model p c = 0.0630965. By comparison of those values with figures ||, || and 
[7| we see clearly that the sample-to-sample fluctuations of the internal energy have a 
maximum at (or at least very near) the multicritical point for each model. It seems 
reasonable to expect similar behaviour in other random models with Z q symmetry. 
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Figure 6. Sample-to-sample fluctuations of the internal energy for the three-state 
Potts model. 
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Figure 7. Samplc-to-samplc fluctuations of the internal energy for the four-state Potts 
model. 



4. Conclusion 

We have studied the ±J Ising model, three-state Potts model and four-state Potts 
model by the numerical transfer matrix method. We have first shown that the sample- 
to-sample fluctuations of the internal energy are, on the Nishimori line, related to the 
specific heat and are expected to have a peculiar behaviour at the multicritical point. 
It was found indeed from numerical calculations that the sample-to-sample fluctuations 
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become anomalously large at or very near the multicritical point in all the models 
we investigated. Since these fluctuations are related to the specific heat by equation 
(8), they should have a singularity at the multicritical point (as the specific heat is 
singular there) although it is difficult to see this effect within our numerical precision. 
The existence of a peak at the multicritical point in conjunction with the existence 
of a singularity, anyway, suggests that the bond distribution and the distribution of 
frustration go under an anomalous change of the behaviour at p — p c . This point 
of view is in agreement with the argument in reference || that the phase transition 
at the multicritical point along the Nishimori line is induced by geometrical anomaly 
H Since the rate of the thermal excitations (e~ /3J /e /3J ) coincides with the ratio of the 
bond configurations ((1 —p)/p) on the Nishimori line, the physical situations are very 
special there, which may be one of the reasons for the existence of a maximum in the 
configurational variance of the thermal quantity at the multicritical point along the 
Nishimori line. It is an interesting future problem to investigate what happens on the 
sample-to-sample fluctuations of physical quantities away from the Nishimori line. 

Note added in proof. We have found after submission of the paper that the sample- 
to-sample fluctuations of the energy, equation (8), is proportional to the derivative of 
the energy by K p with K fixed. Thus a maximum in these fluctuations means that 
the energy changes most dramatically at the multicritical point as K p changes. The 
exact upper bound on the specific heat is the sum of two derivatives by K and K p and 
corresponds to the derivative along the Nishimori line. 
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percolation transition. 



